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Active stabilization in systems with zero or negative stiffness is an essential element of a wide 
variety of biological processes. We study a prototypical example of this phenomenon at a micro¬ 
scale and show how active rigidity, interpreted as a formation of a pseudo-well in the effective energy 
landscape, can be generated in an overdamped ratchet-type stochastic system. We link the transition 
from negative to positive rigidity with correlations in the noise and show that subtle differences in 
out-of-equilibrium driving may compromise the emergence of a pseudo-well. 
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The response of a living system to mechanical loading 
depends not only on the properties of the constituents 
and their connectivity, but also on the presence of non- 
thermal endogenous driving [T]. Thus, molecular motors 
can either stiffen the cytoskeleton through actively gen¬ 
erated pre-stress [2] or fluidize it by facilitating remodel¬ 
ing 0. Powered by ATP hydrolysis, living systems can 
also operate in mechanical regimes with negative passive 
stiffness as in the case of hair cells mu and muscle half¬ 
sarcomeres 00. In those cases metabolic resources are 
used to modify the mechanical susceptibility of the sys¬ 
tem and stabilize the apparently unstable states EHin]. 

At the structural level, active rigidity may be the out¬ 
come of tensegrity tightening m, connectivity change 
[12] , steric interactions m, or the prestress exploit¬ 
ing strong nonlinearity of the passive response mnsj. 
ATP induced stiffening can even take place at the level 
of individual structural elements as in the case of the 
Frank-Starling effect in cardiac muscles that cannot be 
explained by a simple filament overlap change m- 

In this Letter we show that active rigidity can also 
emerge at the micro-scale level through resonant non- 
thermal excitation of molecular degrees of freedom as in 
the case of an inverted pendulum m- Following this in¬ 
ertial prototype, we construct an example of a mechan¬ 
ically unstable overdamped system where stabilization 
and creation of a new pseudo-well in the effective en¬ 
ergy landscape can be induced by a colored noise. The 
proposed mechanism of rigidity generation requires a fi¬ 
nite distance from equilibrium and is therefore different 
from the more conventional entropic stabilization [TS] . 
The possibility of actively tunable rigidity opens inter¬ 
esting prospects not only in biomechanics [T9| but also 
in engineering design incorporating negative stiffness m 
or aiming at synthetic materials stabilized dynamically 

[aiEg. 

We illustrate our idea on a simple bi-stable mechanical 


system described by a single collective variable: the neg¬ 
ative stiffness is viewed as a result of coarse-graining in a 
microscopic system with domineering long range interac¬ 
tions [23|. We assume that this ’snap-spring’ is exposed 
to both thermal and correlated noises and acts against a 
linear spring which qualifies it as a molecular motor oper¬ 
ating in stall conditions [24] . Instead of the conventional 
focus on active force, we study in this Letter a possibility 
of generating by this motor active susceptibility. 

The advantage of our analytically transparent setting 
is that we can distinguish the separate effects of thermal 
(scaled with temperature D) and non-thermal (scaled 
with affinity A) components of the noise on the effective 
energy landscape. We construct a non-equilibrium phase 
diagram in the space of parameters {D, A) showing that 
bifurcations connecting different dynamic phases may be 
either sub-critical, indicating a first order phase transi¬ 
tion, or super-critical, indicating a second order phase 
transition, with the overall behavior controlled by a tri- 
critical point. Some features of the observed dynamic 
transitions are reminiscent of the behavior of the Ising 
model in periodic magnetic field [25] and the behavior 
of folded proteins subjected to periodic forces [26]. We 
show however, that our system is highly sensitive to the 
stochastic nature of the nonequilibrium reservoir. Thus, 
in case of a periodic or dichotomous (DC) rocking force, 
a pseudo energy well exists in an extended domain of the 
parameter space, while it completely disappears if the 
noise is of Ornshtein-Uhlenbeck (OU) type. 

Model. Consider the non-dimensional Langevin equa¬ 
tion X = —dxEP\/2D^{t), where ^(t) is a standard delta 
correlated noise with zero average, D is a measure of tem¬ 
perature and E{x,z,t) = V{x) P k{x — /2 — xf{t) is 

a time dependent potential. We assume that V{x) = 
{\x\ — 1/2) ^/2 is a bi-stable potential describing the con¬ 
formational change and z is a control parameter coupled 
through a spring of stiffness k with the internal variable 
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FIG. 1. Force-elongation curves in the case of periodic driving 
(adiabatic limit). The equilibrium system (A = 0) is shown 
in (a) and and out-of-equilibrium system ^ 0) - in (b). 
Here k = 0.6. 


X. The energy is supplied to the system by a rocking 
force f{t) characterized by an amplitude A and a time 
scale r. The presence of a thermal noise in this mean 
field description suggests that the system has a finite 
size m- The external force required to maintain the 
system in the steady regime is T{z) = k[z — (x)], where 
(x) = lim (1/t) xp(x,t)dxdt is the double averag- 

t^OO j oo 

ing over ensemble and time and p(x, t) is the correspond¬ 
ing probability distribution. 

The main object of our study is the effective potential 
F{z) = f^T(s)ds. While 2 ; has been introduced as a 
constant parameter, it can be also viewed as a mesoscopic 
variable satisfying nz = —dzE + fe^t + V^2z7D^(t) where 
the frictional timescale vjk ^ r and is a slowly 
varying external force. As we show in [28] the dynamics 
of the ensemble and time averaged 2 ; that we denote by Z 
is governed by vZ = —dzF^fext- In the case of skeletal 
muscles, if x characterizes the state of a generic cross 
bridge, Z would be a measure of strain at the level of 
the whole half-sarcomere [28] . Even in the absence of an 
explicit mesoscopic variables, the concept of an effective 
potential is useful for the study of the slow component of 
the motion [29] . 

Periodic driving. Suppose first that the driving is 
square shaped f{t) = with n{t) = [2t/rJ, 

where brackets denote the integer part. An analyti¬ 
cally transparent case is when the correlation time r is 
much larger than the escape time for the bi-stable po¬ 
tential E, e.g. [30]. For f{t) = A we can find the 
stationary solution of the Fokker-Planck equation (with 
2 ; = const) dtP = dx [pdxE + DdxP] analytically, com¬ 
pute {x{A)) explicitly and then average the result over 
the period {{x)) = [{x{A)) + {x{—A))]/2, see [28] for de¬ 
tails. The typical force-elongation curves T{z) and the 
corresponding potentials F{z), obtained in such adia¬ 
batic limit, are shown in Fig. The equilibrium sys¬ 
tem with A = 0 exhibits negative stiffness at 2 ; = 0 
where the effective potential F{z) has a maximum (spin- 
odal state). As temperature increases at A = 0 we ob¬ 
serve a standard entropic stabilization of the configura¬ 
tion 2 : = 0, see Fig. [^a), which takes place as a second 




FIG. 2. Phase diagram in (A, D) plane showing phases I,II 
and III: (a) - adiabatic limit, (b) - numerical solution at r = 
100 (b). Ca is the tri-critical point, De is the point of second 
order phase transion in a passive system. The Maxwell line 
for first order phase transition is shown by dots. Here k = 0.6. 


order phase transition at the equilibrium temperature 
De = r/[8(l + k)] where r is a root of a transcenden¬ 
tal equation 1 + yJVJi^e~^Fj\^ ^ erf(l/v^)] = rf{2k) 
[28] . As the degree of non-equilibrium, characterized by 
A, increases, the effective potential develops a pseudo¬ 
well with a minimum at 2 : = 0, see Fig. [^b), and we 
associate this phenomenon with the emergence of active 
rigidity. 

The non-equilibrium steady state (dynamic) phase di¬ 
agram summarizing the results obtained in adiabatic ap¬ 
proximation is shown in Fig. (a). There, the ’param¬ 
agnetic’ phase I describes the regimes where the effective 
potential F{z) is convex, the ’ferromagnetic’ phase II is a 
bi-stability domain where the potential F{z) has a double 
well structure and, finally, phase III is where the function 
F{z) has three convex branches separated by two concave 
(spinodal) regions. If we interpret the boundary Ca — De 
separating phases I and II as a line of (zero force) second 
order phase transitions and the dashed line Ca — Ma 
as a Maxwell line for the (zero force) first order phase 
transition, see [28], then Ca will be a tri-critical point. 
Near this point the system can be described by the non- 
equilibrium Landau potential F{z) = + 

where the coefficients r, q are the measures of passive and 
active excitations, respectively, while p > 0 is a fixed pa¬ 
rameter. Similar tri-critical point has been observed in 
the periodically driven mean field Suzuki-Kubo model of 
magnetism [31] which can be interpreted in our terms as 
a description of the T = 0 behavior only. 

The adiabatic approximation fails at low temperatures 
(small D) where the escape time diverges and in this do¬ 
main the corrected phase diagram was obtained numeri¬ 
cally by computing the appropriate periodic solutions of 
the Fokker-Plank equation, see Fig.[^(b). The high tem¬ 
perature part of the diagram (tri-critical point, point De 
and the vertical asymptote of the boundary separating 
phases I and III at large values of A are captured ade¬ 
quately by the adiabatic approximation. The new feature 
is a dip of the boundary separating Phases II and III at 































3 


some D < De leading to an interesting re-entrant behav¬ 
ior (cf. [32j[33]) which is an effect of stochastic resonance. 
To verify our numerics in the low temperature domain we 
used the Kramers approximation, to show that indeed 
A = 1/2 at point K and 74 = l/2 + /c/4at point M, see 

m- 



FIG. 3. (a-c) Typical tension-length relations in phases I, II 
and III. Points a, (3 and 7 are the same as in Fig. [ 2 ](b); (d) 
shows the active component of the force. Inserts illustrate the 
behavior of stochastic trajectories in each of the phases (gray 
lines) and their ensemble averages (black lines) at ~ 0. The 
other parameters: k = 0 . 6 , r = 100 . 

Force-elongation relations in different points of the 
(A^D) phase diagram (Fig. (b)) are shown in Fig. 
where the insets illustrate the typical stochastic trajec¬ 
tories. We observe that while in phase I thermal ffuctu- 
ations dominate periodic driving and undermine the two 
wells structure of the potential, in phase III the jumps 
between the two energy wells are fully synchronized with 
the rocking force. In phase II the system shows inter¬ 
mediate behavior with uncorrelated jumps between the 
wells. We conclude that the pseudo-well in phase III 
has a resonant nature and remark that somewhat similar 
phenomena were also observed in other driven out-of- 
equilibrium systems [34]. 

In Fig.j^d) we show the active component of the force 
Ta{z) = T{z;A) — T{z;0) representative of phases I, II 
and III. The active contribution is significant only in 
phase III and the corresponding plateau can be viewed as 
another signature of the presence of a pseudo-well. Inter¬ 
estingly, our prototypical device generates active tension 
of both signs which can be interpreted as pulling at z > 0 
and pushing at z < 0. However, in the puling regime 
the linear spring is stretched while in the pushing regime 
it is compressed. Since in biological conditions the fil¬ 
aments responsible for passive stiffness would buckle in 


compression, e.g. [35], the pushing part of the active 
force-length relation is hardly realistic. On the other 
hand, the pulling part shows a striking resemblance to 
the isometric tetanus in skeletal muscles [36] that can be 
also driven through the bi-stable potential [37] . 



FIG. 4. (a) Phase diagram in the case of DG driving. The 
identification of phases I, II,III is the same as in Fig. 11 (a,b). 
(b) Typical tension-length relations in different phases (b). 
Here r = 100 and k = 0.6. 

In view of this analogy, detailed in [28] , it is instructive 
to estimate the four non-dimensional parameters of the 
model by using the available data on molecular motors 
operating in muscle cells. We choose the time scale to 
he = r]/ko 0A ms where 77 ^ 0.38 ms. pN/nm is 
the viscosity adopted in [7] and ko ^ 3 pN/nm is the 
stiffness of the cross-bridge in pre and post power stroke 
configurations. The spatial scale is /* = a ^ lOnm, the 
characteristic size of a motor power-stroke [40] and the 
stress scale is k'^ = ko. This leads to an energy scale 
= koa^ ^ 200 pN.nm. Then, the non-dimensional 
parameters can be estimated as follows. Parameter 
k = kmjk'^ ^ 0 . 6 , where km is the stiffness of the elas¬ 
tic part of the myosin motor [38] |39]. Temperature is 
D = ksO/E'^ ^ 0.01 where ks = 4.10 pN.nm is the 
Boltzmann constant and 0 ^ 300K the ambient tem¬ 
perature. For the active driving time scale, we esti¬ 
mate r = Tajr'^ ^ 100 where = 40 ms is the char¬ 
acteristic time of ATP hydrolysis m- Finally we take 
A = y/Ajj./E'^ ^ 0.5 where A/i = 20kB& is the degree of 
non-equilibrium of the hydrolysis reaction [41]. The ob¬ 
tained estimate {A = 0.5,1) = 0.01) suggests that mus¬ 
cle myosins, operating in stall conditions (isometric con¬ 
tractions), are in phase HI. The proposed representation 
of the ATP hydrolysis (through parameter A) explains 
stabilization of the power stroke mechanism in skeletal 
muscles in the negative stiffness regime [7] and may be 
also behind titin based force generating mechanism at 
long sarcomere lengths that does not rely on actin-myosin 
based cross-bridge interactions [42] . 

Dichotomous driving. To ascertain the robustness of 
these results we now consider a different representation 
of the external forcing as a dichotomous (DC) or tele¬ 
graphic noise, e.g. [l^lll]. In this case f(t) = A(— 1 )^^^\ 
where n{t) is a Poisson process with {P{n) = e~^X^/nl 
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and rate parameter A = l/(2r); we thus have (/(t)) = 
Ae'xp{—t/r) and (/(t),/(s)) = ^4^ exp(—|t — s|/r). The 
DC driven system is controlled by the same number of pa¬ 
rameters as the periodically driven system, however, the 
problem is no longer analytically tractable. The numeri¬ 
cal solution of the ensuing stochastic differential equation 
shows that the qualitative structure of the phase diagram 
in the (^4, D) plane remains the same as in the case of 
periodic driving, see Figll We checked our numerical re¬ 
sults by considering an analytically tractable double limit 
when r ^ 0, ^4 ^ oo, while D = A^r remains finite. In 
this limit phase III disappears because the system can 
be viewed as exposed to a white noise with effective tem¬ 
perature D* = \/ D‘^ D‘^. Then there is only a second 
order phase transition at the expected value of the pa¬ 
rameter D* = r/[ 8 (l + k)]. This simple limit highlights 
the crucial role of correlations in the noise (r 7 ^ 0). Our 
next example, however, shows that correlations per se are 
not enough. 

Ornstein-Uhlenheck driving. Suppose now that f{t) is 
a solution of a linear stochastic differential equation / = 
where Cf{t) is a standard white 
noise independent of ^(t). Such non-equilibrium driving 
is known as Ornstein-Uhlenbeck (OU) noise, e.g. [441145] 

, and its first ((/(t))) and second moments 

are the same as in the case of DC if we assume, without 
loss of generality, that /(O) = A. The Fokker-Planck 
equation for the probability density p(x, /, t) takes the 
form dtP = d^ipd^E + Dd^p) + T~'^df{fp + A^dfp). By 
solving it numerically we obtain a phase diagram shown 
in Fig. I^a). A striking feature of this diagram is that 
phase III is missing because, in contrast to periodic and 
DC case, the noise is now unbounded and the system can 
always escape from a neighborhood of a resonant state. 
The behavior of the force-elongation relations shown in 
Fig. I^b) is compatible with the idea of purely entropic 
stabilization, in particular, the limit of thermal noise is 
again recovered when r —0 and A ^ oc, with D = A^r 
fixed. 



FIG. 6. (a) Phase diagram in the case of DC driving. The 
identification of phases I, II and III is the same as in Fig. 
(a,b). The dash-dotted line - the boundary between phases 
II and III (periodic driving); the dashed line - the boundary 
between Phases I and II (OU driving), (b) Typical tension- 
length relations in the case of DC driving in different phases 
(b). Parameters k = 0.6, D = 0. 


ferences between our three representations of a non¬ 
equilibrium bath, we compare in all three cases the (r, D) 
phase diagrams corresponding to the limit D ^ 0 where 
the thermal component of the noise is absent. In the DC 
case the solution of the limiting Fokker-Plank equation 
can be written explicitly [46] 

Pnc(^) ~ <3(x)-‘ exp 1 ^ 

where Q = A^ - (5^U(x))2 and U = U + k{x - zf/2. 
The choice of the normalization constant depends on the 
parameters and is detailed in [28|. The resulting phase 
diagram, shown in Fig. [^a), exhibits all three phases 
with a tri-critical point C located at tc = [ 2 (/c + 1 )]“^ 
and Dc' = De + [2{k + l)]“^/4. The behavior of the 
force-elongation relations in different phases is illustrated 
in Fig. [^b). 

In the case of OU driving with D = 0 an analytical 
approximation of the stationary probability distribution 
is available for small r only [46] 




FIG. 5. (a) Phase diagram in the case of OU driving. The 
identification of phases I, II is the same as in Fig. [2](a,b). (b) 
The typical tension-length relations in different phases. Here 
T = 100 and k = 0 .6. 


Zero temperature limit. To clarify further the dif¬ 


/ ^ ^ ( V(x)^Tid^V(x)Y/2\ 

Pou{x) - R{x) exp I- V U / j 

where R{x) = \l-\-Tdxxy {x)\^ see [28] for the details. The 
resulting phase diagram does not contain phase III and 
the line dividing phases I and II is shown in Fig. [^a) 
(dashed line). The problem with periodic driving ex¬ 
hibits in the limit D —> 0 only phases II and III even for 
rapidly oscillating external fields, see the dash-dotted line 
in Fig.j^a). In this perspective, the DC driving emerges 
as an intricate amalgam of OU and periodic noises with 
none of them dominating the other. 

Conelusions. To complement the existing microscopic 
models of force generation (Brownian ratchets), we pro¬ 
posed a conceptually similar model of rigidity genera¬ 
tion (Brownian snap-springs). The model, invoking some 
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interesting parallels between condensed matter physics 
and biomechanics, shows that by controlling the degree 
of non-equilibrium in the system, one can modify the 
structure of the effective energy landscape. In particular, 
this implies that unstable or marginally stable mechan¬ 
ical states may be stabilized by out-of-equilibrium ATP 
hydrolysis reaction. Our results also suggest that the 
mechanical action of a non-equilibrium reservoir can be 
crucially sensitive to the higher moments of the stochas¬ 
tic noise. 

The authors thank J.-F. Joanny, R. Garcia Garcia and 
M. Camel for helpful discussions. 
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